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Abstract. We have studied the evolution with temperature of the local as well as the average crystal 
structure of MgB2 using the real-space atomic pair distribution function (PDF) measured by high resolution 
neutron powder diffraction. We have investigated the correlations of the B-B and B-Mg nearest neighbor 
pair motion by comparing, in the wide temperature range from T — 10 K up to T — 600 K, the mean- 
square displacements (MSD) of single atoms with the mean-square relative displacements (MSRD) obtained 
from the PDF peak linewidths. The results show that the single atom B and Mg vibrations are mostly 
decoupled from each other, with a small predominance of positive (in phase) correlation factor for both 
the B-B and B-Mg pairs. The small positive correlation is almost temperature independent, in contrast 
with our theoretical calculations; this can be a direct consequence of the strong decay processes of the E2 g 
anharmonic phonons. 

PACS. 74.70. Ad Metals; alloys and binary compounds (including A15, MgB2, etc.) - 61.T2.-q Neutron 
diffraction and scattering - 63.20.-e Phonons in crystal lattices 

pling is mainly driven by the interaction between elec- 
tronic carriers in the 2D a band with boron p XtV character 
and the zone center E 2g phonon mode. [T^ll3ll4ll5ll6ll7| 
This is reflected in a Kohn anomaly in the phonon dis- 
persion related to the size of the small 2D tubular Fermi 
surfaces. The proximity of the Fermi level to the Van 
Hove singularity (VHs) and to the band edge discloses a 
new scenario where the large amplitude of the expected 
boron zero point lattice fluctuations |7I12I15I18I21I22| in- 
duces large fluctuations of the same order of the separation 
between the VHs, the gap edge and the Fermi level itself. 
Although the amplitude of the lattice fluctuations seems 
thus to be highly relevant for the superconductivity in 
MgB2 there is a lack of experimental information on this 
key point. Furthermore, although the average structure 
(P6 /mmm) of the MgB 2 system has been exhaustively in- 
vestigated, there is not yet any study of the local structure, 
since typical x-ray local probes, such as EXAFS, cannot 
be used to study local structure near light atoms. 

In view of this, here we employ high resolution neutron 
diffraction to obtain the pair distribution function (PDF) 
of MgB2- In this way we investigate the local as well as the 
average structure of MgB2 , namely the mean-square dis- 
placements (MSD) of single atoms and the mean-square 
relative displacements (MSRD). The comparison of these 



1 Introduction 

MgB 2 is the simplest system to investigate the quantum 
mechanism that allows the formation of a superconduct- 
ing condensate with critical temperature T c ~ 40 K pQ a 
factor two higher than in all other known intermetallic su- 

\ perconductors. It was recently proposed by a few groups 
[2I4I3I5I5] that the enhancement of the critical tempera- 
ture in MgB 2 is due to the exchange-like interband pairing 

, in a multiband superconductor. There is now experimen- 

■ tal evidence |7l8l9ll(JflT] that MgB 2 is the first clear case 
of a high T c multiband superconductor showing two-gaps 
in the a and n bands respectively, in agreement with the 

' theory. [HSI 

The characteristic feature of MgB 2 is that the electron- 
phonon interaction gives a weak pairing in the 7r channel 
and a strong pairing in the cr-channel. Electronic band 
calculations. [141151161171 Raman [16119) and inelastic 
neutron scattering experiments |2l)| provide evidence of 
an extremely large deformation potential for the B bond 
stretching modes, which gives rise to strongly anharmonic 
phonons. This anharmonicity also results in a structural 
instability (phonon softening) that affects the dynamics 
of the lattice fluctuations and the local structural proper- 
ties of the material, as we discuss below. There is now a 
general agreement that this strong electron-phonon cou- 
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quantities permits for the first time to extract the corre- 
lation factors pb-b, PB-Mgj defined below, of the boron- 
boron and boron-magnesium pair motions, which are found 
to be pb-b ~ 0.1, pB-Mg ~ 0.1, and nearly constant 
in a wide range of temperature K < T < 600 K. We 
also compare the experimental data with a constant force 
(CF) model for the phonon dispersion. We estimate that 
the phonon frequency renormalization due to the electron- 
phonon interaction on the E^ g modes yields a reduction 
ApB-B ~ —0.03 in the boron-boron correlation factor. 
While the CF model can nicely account for the zero tem- 
perature values of the single atoms MSD and the corre- 
lated MSRD, the temperature behavior of the correlation 
factor is shown to be highly anomalous and its physical 
interpretation gives rise to new questions about our un- 
derstanding to the local lattice dynamics in this material. 



2 Experimental method and data processing 

Polycrystalline samples of MgB2 were synthesized at high 
temperature by direct reaction of the elements in a tan- 
talum crucible under argon atmosphere using pure n B 
isotopes. Time-of-flight neutron powder diffraction data 
were collected on the NPDF diffractometer at the Manuel 
Lujan, jr., Neutron Scattering Center (LANSCE) at Los 
Alamos National Laboratory. The NPDF diffractometer, 
with its high neutron flux and backscattering detector 
modules has a high resolution, Ell providing access to a 
wide range of momentum transfer with sufficient counting 
statistics making it ideal for PDF studies. The powdered 
samples were sealed inside extruded cylindrical vanadium 
containers. These were mounted on the stage of a cry- 
ofurnace with and without heat-shield for T < 300 K 
and T > 300 K, respectively. The scattering data from 
the empty cryofurnace, with and without heat-shield, an 
empty container mounted on the cryofurnace, and the 
empty instrument were also collected, allowing us to as- 
sess and subtract instrumental backgrounds. The scatter- 
ing from a vanadium rod was also measured to allow the 
data to be normalized for the incident spectrum and de- 
tector efficiencies. We collected each diffraction spectrum 
up to the high momentum transfer of Q = 40 A^ 1 in 3 
hours. The high resolution diffraction spectrum so ob- 
tained presents both the Bragg peaks and the diffuse scat- 
tering. While the Bragg peaks reflect the long-range or- 
der of the crystalline samples, the oscillating diffuse scat- 
tering contains local structural information including the 
correlated dynamics. |29l31j The PDF is obtained from 
a Fourier transform of the powder diffraction spectrum 
(Bragg peaks + diffuse scattering). [23 It consists of a se- 
ries of peaks, the positions of which give the distances of 
atom pairs in the real space, while the ideal width of these 
peaks is due both to relative thermal atomic motion and 
to static disorder. This permits the study of the effects 
of lattice fluctuations on PDF peak widths and yields in- 
formation on both the single atom mean-square displace- 
ments and the relative mean-square displacements of atom 
pairs and their correlations. 



Standard data corrections |23| were carried out using 
the program PDFGETN.E5| After being corrected, the 
data were normalized by the total scattering cross sec- 
tion of the sample to yield the total scattering structure 
function S(Q). Afterwards, the total scattering structure 
function S(Q) is converted to the PDF, G(r), by means 
of a sine Fourier transform according to the relation: 

G(r) = 47rr( ( o(r) - p ) 

= ~ / Q[S{Q)-l]sm{Qr)dQ. (1) 
n Jo 

We modeled the PDF using a structural model that takes 
advantage of the definition of the radial distribution func- 
tion RDF R(r), namely: 

R(r) = rG(r) + 47rr 2 p = £ &{r r«), (2) 

where bi is the scattering length of the i th atom, (b) is the 
scattering length averaged over the sample composition, 
rij = |r, — Tj\ is the distance between the i th and the 
j th atoms, and the sums are taken over all the atoms in 
the sample. Before being compared to the data, the cal- 
culated G(r) is convoluted with a termination function, 
sm(Q max r)/V to account for the effects of the finite data 
collection range. Fundamental lattice information, such as 
the average crystal structure, the lattice constants, a scale 
factor, and the refined atomic (thermal) displacement pa- 
rameters, can now be extracted from the PDF by using 
the PDF refinement program PDFFIT E5] that is based 
on a least-squares approach to fit the PDF profile. The av- 
erage atomic displacement distribution of atom i along the 
major axes x, y, and z, a 2 (i x ), a 2 (i v ), a 2 (i z ), are defined 
as 

a 2 (i a ) = ([ Ul -f Q ] 2 ), (3) 

where is the lattice displacement of atom i from its 
average position j7| and r a is the unit vector pointing 
along the direction a = x,y, z. Due to the geometry of 
these compounds, the two boron atoms for the unit cell 
are equivalent and a 2 (i x ) = <J 2 (i y ) for all the atoms, so 
that only four parameters were needed namely a 2 (B xy ), 
a 2 (B z ), a 2 (Mg xy ) and a 2 (Mg 2 ). 



3 Results and discussion 

In Fig.^we show the reduced scattering structure function 
Q[S{Q) - 1] for the MgB 2 at T = 300 K, while the cor- 
responding reduced PDF, G(r), obtained using Eq. (JIJ, is 
shown in Fig. [21 The features of the NPDF diffractometer 
allowed us to obtain high quality PDFs as can be noted 
by inspecting the modeled fit (solid line) of the G(r) in 
the upper panel Fig. E] This can be seen also in the lower 
panel of Fig. EJ where most of the fluctuations in the differ- 
ence curve (solid line) are within the standard deviation 
on the data ±Z\[G(r)] (dashed lines above and below the 
difference curve). 
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frequency dispersions are expected to be notably differ- 
ent from dispersionless Einstein models in these materials. 
Moreover, the Einstein fits were obtained by simply using 
an effective atomic mass Mm s = 24.3 a.m.u. for a 2 (Mg) 
and Mb = 10.81 a.m.u. for er 2 (B). This means that we 
are implicitly assuming that magnesium and boron vibra- 
tions are totally decoupled from each other. These two- 
parameter fits reproduce the data very well for each of 
the motions considered. Non-independent boron and mag- 
nesium vibrations would lead to deviations from Einstein 
behavior that are not evident in the data suggesting that 
this approximation is reasonable in this system. We return 
to this point below. 

To gain further insight on this issue we introduce a 
constant force (CF) shell model for the dynamical matrix. 
We neglect for the moment the effects of the electron- 
phonon interaction, which leads to rcnormalizcd phonon 
frequencies and anharmonicity, and we assume the lattice 



3.1 Atomic mean-square displacements 

In order to extract information about the phonon spec- 
trum, we fit each lattice displacement with a simple Ein- 
stein model, 



a 2 (i a ) 



MiLOi, 



+ Co(*a), 



(4) 



where 0o(i Q ) takes into account the static disorder, Mi 
is atomic mass of the atom i and tu(i a ) provided an es- 
timate of the vibrational frequency of the atom i along 
the direction a. The quantity n(x) is the Bose thermal 
factor n[x] = [exp(x/k B T) — 1] . The values of the fit- 
ting parameters, Wi a and cto(j q ), are reported in Tab. ^ 
and the a 2 vs. T fitting curves are represented by dotted 
lines in Fig. [3J The different values of Lo ia represent the 
different energy range of the phonon spectra associated 
with the boron and magnesium in-plane and out-of-plane 
lattice vibrations, and they are in good agreement with 
the corresponding spectra reported in Ref. 

Quite surprising, we find that, apart from some amount 
of disorder in cr 2 (Mg z ), this simple four-peak Einstein 
model seems to describe quite well the lattice vibrations 
in this compound. This is quite intriguing because phonon 



The PDF's have been fit over the range 1 — 18 A using a 
hexagonal crystal structure (space group P6/mmm), and 
excellent fits were obtained at all temperatures. The Bragg 
peaks are clearly persistent up to 25 A -1 , reflecting both 
the long-range order of the crystalline samples and the 
small amount of positional (vibrational or static) disorder 
of the atoms around their average positions. 

The lattice parameters at T = 300 K for MgB 2 (T c ~ 
39 K) were found to be a = 3.08505(4) A, c = 3.5218(1) A. 
The lattice displacements a 2 (B xy ), <r 2 (B z ), a 2 (Mg xy ), cr 2 (Mg z ) 
are shown in Fig. [3] (open circles) . 
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Fig. 2. (Upper panel) The PDF G(r) = 4?rr[p(r) - p ], ob- 
tained from Eq. 0for MgB2, measured at 300 K (dots), with 
the structure refinement curve obtained by a least-squares ap- 
proach (solid line). (Lower panel) The difference curve (solid 
line) of the experimental PDF with the modeled fit and the 
standard deviation on the data ±Zi[G(r)] (dashed lines) are 
shown. We can observe that most of the fluctuations in the 
difference curve are within ±zi[G(r)]. 
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Fig. 3. Anisotropic in plane a 2 (B xy ), a 2 (Mg xy ) and along 
the c-axis <t 2 (B z ), u 2 (Mg z ) mean-square displacements for the 
B and Mg atoms, as a function of the temperature. The open 
circles represent the PDFFIT refined values, while the solid 
and the dotted lines represent the modeled data in the CF and 
Einstein model, respectively. The errorbars in the a 2 (i a ) are 
smaller than the size of the used symbols (open circles). The 
dashed line shows the in-plane boron lattice displacements in 
the CF model in the absence of the electron-phonon interac- 
tion. 



u{B xy ) (K) 


703 ± 10 


w(B„) (K) 


572 ± 15 


w(Mg*J (K) 


398 ±6 


w(Mg.) (K) 


350 ±9 






^(B»„) (A a ) 


< 5 ■ 10~ 5 


°o(B z ) (A 2 ) 


< 9 ■ 10~ 5 


*S(M|U (A 2 ) 


< 7 ■ 10" 5 




0.0013 ± 0.0003 



Table 1. "Effective" phonon frequencies uJi a and static dis- 
order contributions a 2 (i a ) for the in-plane and out-of-plane B 
and Mg displacements as obtained by the Einstein model fit 
[Eq. 0] of a 2 (i a ). 



dynamics to be harmonic (we shall discuss later and more 
specifically the role of the electron-phonon interaction and 
possible effects of the anharmonicity of the Ei g phonon 
mode). This will enable us to evaluate eigenvectors e qiAt 
and eigenvalues w„„ of the lattice modes for each point of 
the phonon Brillouin zone. The phonon contribution to the 



lattice displacements a 2 (i a ) 
will be thus obtained as:j!T3] 



ph 



(ia) + <?o(ia)} 



fph(*a) 



N 2-*> 



MiU! n 



- + n(cj q)ft ) 



(5) 



4>r (ev/A j 


H (eV/A ) 


<p± (ev/A ) 


12,45 


49.80 


21.17 


Xr (eV/A 2 ) 


XII (eV/A 2 ) 


XX (eV/A 2 ) 


0.0 


16.6 


0.0 


K r (eV/A 2 ) 


«H (eV/A a ) 


k ± (eV/A^) 


0.84 


1.50 


9.14 


i> r (eV/A a ) 


^11 (eV/A a ) 


^x (eV/A a ) 


0.0 


9.34 


0.42 



Table 2. Force constant parameters reproducing the bare 
phonon dispersion in MgB2 in the absence of electron-phonon 
interaction. 



where is the component of the eigenvector e qiA1 con- 
cerning to the displacement of the i atom along the a 
direction and N is the total number of q-points consid- 
ered in the phonon Brillouin zone. From a general point 
of view, since cr 2 h (i a ) involves an integral over the whole 
Brillouin zone and over all the phonon branches, it will not 
be sensitive to the fine details of the phonon dispersion but 
only to its gross features. For this reason, and in order 
to preserve the simplicity of our analysis, we limit our- 
selves to consider only four elastic springs, 4>, x, k and ip, 
connecting, respectively, in-plane B-B nearest neighbors, 
out of plane B-B nearest neighbors, out-of-plane B-Mg 
nearest neighbors, and in-plane Mg-Mg nearest neighbors. 
Each elastic spring is specified by its tensor components 
(ex.: <j) r , </>||, (/>j_), corresponding respectively to the lattice 
displacements along the radial (bond-stretching) direction 
and along the in-plane and out-of-plane tangential (bond- 
bending) directions. In MgB2 we choose the constants <fi r , 

0|| , <P±, Xr, X\\, X±, *r, K h K -L' Vv, ip±, to fit the 
local-density functional (LDA) phonon dispersion of Ref. 
Hi along the high-symmetry points of the Brillouin zone. 
Since the constant force model is meant to reproduce the 
bare phonon dispersion, we deliberately did not include in 
the fit procedure the E2 g phonon frequencies at the zone 
center r, A, which are known to be strongly affected by 
the el-ph interaction. The elastic constants obtained from 
this fitting procedure are reported in Tabled 

As mentioned in the introduction, the Ei g phonon 
modes close to the F and A points are expected to be 
strongly affected by the interaction with the almost 2D 
parabolic a bands, giving rise to a remarkable softening of 
the E 2g phonon frequencies for |q| < 2fcp, where fcp is the 
Fermi vector of the a bands. [T4T15I16I17I2(J| We include 
these effects through the self-energy renormalization of the 
phonon frequencies i7| 2S (q) = w| 2g (q)-(4A^ <T g 2 / anharm /M B )il2D(q), 
where N„ is density of states of the a bands per spin and 
per band, g the electron-phonon matrix element between 
cr-band electrons and the Ei g phonon mode at the zone 
center and /anharm is a dimensionless factor accounting 
for the anharmonic hardening of the Ei g phonon modes 
due to the electron-phonon coupling itself. Moreover the 
factor 4 takes into account the spin and band degener- 
acy and -/72D(q) is the two-dimensional Lindhardt func- 
tion 77 2D (x) = 0{l-x)+9(x~l)[x-y / x 2 - l]/a; 4 , with x = 
|q|/2fcfp. We take, from first-principle calculations. |14ll5H8) 
/anharm = 1-25, N a = 0.075 states/(eV • spin • cell) g = 12 
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Fig. 4. Phonon dispersion and phonon density of states for 
MgB2 evaluated within the CF model. Dashed lines represent 
the same quantities without taking into account the electron- 
phonon interaction. 



eV/A and hp ~ tt/12cIb~b, where de-B is the boron- 
boron distance. 

The phonon dispersion and phonon density of states 
(PDOS) of our constant force model are shown in Fig. 
["fl in good agreement with the LDA calculations of Ref. 
|lfij . For comparison, the dashed lines represent the bare 
phonon dispersion in the absence of el-ph interaction and 
the corresponding PDOS. Most striking is a partial shift 
in spectral weight from 100 meV to 70 meV when the el- 
ph coupling is turned on, due to the softening of the flat 
r — A band originating from the i?2 ff phonons. 

The lattice displacements for each atom and for each 
direction can now be evaluated directly by means of Eq. (J5J • 
The results are also shown in Fig. leaking with the data ex- 
tracted by the PDF analysis and with the simple Einstein 
fits. In spite of the crude simplification of the phonon dis- 
persion, the agreement with the PDF data is remarkable. 
Just as in the Einstein fits, we find that a constant shift 
ffp(Mg z ) ~ 0.0013 A 2 due to intrinsic disorder is needed to 
account for the magnesium out-of-plane lattice vibration 
while the contribution of the disorder for the other modes 
is found to be negligible. Note that the agreement between 
the data and solid lines was obtained with no adjustable 
parameters. The values used in the constant force mod- 
els in Table |21 were indeed fit to reproduce the theoretical 
LDA calculation and not the PDF data and the only other 
parameter, the static disorder parameter in er 2 (Mg z ), was 
the same value as was used in the Einstein fits. The agree- 
ment with the data and the Einstein model fits is also 
very good suggesting that boron and magnesium lattice 
vibrations, as well as in-plane and out-of-plane lattice vi- 
brations, are on average independent each other. 

To point out the explicit role of the electron-phonon 
interaction, we show also in the left-upper panel of Fig. 
the a^(B xy ) lattice displacements in the absence of the 
el-ph frequency renormalization. Since only Ei g phonons 
are coupled, only the <Jq(B xv ) lattice displacements result 
modified. We note that the inclusion of electron-phonon 
interaction effects leads to a slight increase of the amount 



of the boron in-plane lattice displacements, with a better 
agreement with the experimental data. The increase of 
Co (Bxy) is easily understandable as due to the softening 
of the E2 g phonon mode. On the other hand, since the 
electron-phonon renormalization effects are restricted to a 

small region ^Jq%. + < 2fcp of the whole Brillouin zone, 

the impact of the el-ph coupling on the total amount of the 
lattice displacements a^{Q X y) is relatively weak. As we are 
going to see, the effects of the electron-phonon interaction 
are more apparent in the correlated pair motion. 



3.2 Correlations in the B-B and B-Mg atomic pairs 
motion 

Above we showed that the average uncorrelated thermal 
motions (equivalent to the Debye- Waller factor in crys- 
tallography) measured from MgB 2 are well explained by 
harmonic models with independent boron and magnesium 
motions. A strength of the PDF technique is that it is sen- 
sitive to correlations in the atomic dynamics that contain 
some additional details about the underlying interatomic 
potentials. 27 29 30 31 Here we explore the motional cor- 
relations in the MgB2 PDF data. 

The Gaussian width er,j of the PDF peaks is directly 
related to the mean-square relative displacement of atomic 
pairs projected onto the vector joining the atom pairs. |2*7] 
Explicitly, 

~2 _ /ry„ „ \ t. l2\ 



,0 



(G) 



where and Uj are the lattice displacements of atoms i 
and j from their average positions, f ^ is the unit vector 
connecting atoms i and j, and where the angular brackets 
indicate an ensemble average. [2*7) 

In our analysis, we focus on the width of the near- 
est neighbor boron-boron PDF peak, (cr B _ B ), and on the 
nearest neighbor magnesium-boron peak, (c B _ Mg ). These 
are well resolved single-component peaks in the PDF whose 
width directly yields correlated dynamical information. 

cr B _ M , as measured by the 



The Gaussian widths c B _ B , 



PDF data are shown in the upper panels of Fig. U3 (open 
circles) along with the same quantity evaluated within the 
CF model (solid line) according to the relation: 



N 



E 



i 



2Re[(e^ I -f. iJ )(e^-f u )e^ 



(7) 



In this latter case in the right panel for c B _ Mg we have 
also added a small contribution of the local lattice dis- 
placements due to the disorder (see below for more de- 
tails). 

Here we first note that, while the zero temperature 
values of <J B _ B , cr|-Mg are wen reproduced by the CF 
model, this model is less good than was the case for the 
uncorrelated motions, especially at high temperatures. 
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Fig. 5. Mean-square relative lattice displacements al-B ( a )> 
°B-M g (b), as extracted from the width of the PDF peaks 
and the corresponding correlation function pb-b (c), pB-Mg 
(d) (open circles). Solid lines represent the same quantities 
evaluated within the CF model, while dotted lines in the upper 
panels show the predicted relative lattice displacements cr|_ B , 
°B-M g f° r completely uncorrelated lattice motion pb-b = 0, 
PB-Mg = 0. The dashed lines in the upper panels represent 
the relative lattice displacement o" B _ B an< i ^B-Mg evaluated 
within the CF model in the absence of the electron-phonon 
interaction. 



The temperature dependence of the single-atom mo- 
tions was well explained by the model, but not those of the 
correlated <7y, which suggests that the model is not cap- 
turing some aspect of the motional correlations. Introduc- 
ing the electron-phonon coupling into the model improves 
the agreement slightly, with a larger effect observed on 
°b-b> but this does not explain all of the discrepancy. In 
each case, the pair correlation peaks broaden more quickly 
in the data than in the model. We note that the models 
are harmonic. Even in the case where we have introduced 
electron-phonon coupling, harmonic spring constants have 
been obtained from a fit to the LDA bands and we specu- 
late that the discrepancy at high temperature is is a result 
of anharmonicity in the boron motion. 

As a general consideration, we would like to stress once 
more that while o~ 2 {i a ) probes the absolute magnitude of 
the single atom mean-square displacement, afj provides 
information about the correlation between the lattice dis- 
placements of atom pairs. Let us consider for instance 
the case of c B _ B which involves only boron in-plane lat- 
tice fluctuations. We can identify three limiting behav- 
iors for this quantity: i) perfectly in-phase lattice mo- 
tion: ii) perfectly opposite-phase motion; Hi) completely 
independent motion. In the first case it is easy to see 
that c B _ B — 0, while <t b b — 2o- 2 (B xy ) when the near- 
est neighbor boron lattice displacements are uncorrelated, 
Aa 2 (B xy ) when they have opposite phase. 



[(Ui 



[(U r ?y] 2 )-2((U r ?i 



•(Uj-fy 



)>• (8) 



Here the first two terms are related to mean-square ther- 
mal displacement of atoms i and j projected along f^, 
while the third term is a displacement correlation function, 
which carries information about the motional correlations. 
It is now useful to quantify the degree of correlation by 
introducing the dimensionless correlation parameter pij 
defined as:ESESEHS3 



+ <7 2 {ji) - 2a(ij)a(ji)pi_ 



(9) 



where cr 2 (ij) — ([(u^ • ry] 2 ). Positive values of p > de- 
scribe a situation where the atoms move in phase, so that 
the resulting value of a 2 j is smaller than for the uncorre- 
lated case. On the other hand, a predominance of opposite 
phase atomic vibrations should result in p < and in a 
PDF peak width of- larger than the uncorrelated case. 

It is important to note that the correlation function pij 
in Eq. 10 expresses the degree of correlation between the 
total atomic displacements. In the presence of two different 
sources of lattice displacements (phonons and disorder), it 
is more convenient to split a 2 j into a phonon and a disorder 
contribution. Assuming the local lattice displacements due 
to the disorder to be uncorrelated, we can write thus: 



-o- 2 (ij) + a 2 (ji), 



(10) 



where pij represents now only the correlation between 
phononic lattice displacements. 

Using Eq. I|10|l the correlation parameter can be cal- 
culated from the total width of the PDF peak as 



Pi] 



r ph(*i) + CT ph0i) + + a 2 {ji) - a. 



2cr p hfe)o- p h(ji) 



(11) 



and erj|_ B 



Finally, the projected atomic mean-square displacements 
o~ 2 (ij) (lattice vibrations along the pair fy- direction) can 
be related to <t 2 q (lattice vibrations along the Cartesian 
axes) by simple geometrical considerations. We have thus 
^ph(B B ) = ^ h (B xy ), a 2 ph (B Mg ) = [AR 2 a 2 ph (B xy )+3a 2 ph (B z )]/(m 2 
3), a 2 ph (Mg B ) = [AR 2 a 2 ph (Mg xy ) + 3a 2 ph (Mg z )}/ (4R 2 + 3), 
where R = a/c = 0.88 and where a and c are the in-plane 
and out-of-plane lattice constants. Similar relations hold 
true for the disorder contributions. 

The phonon correlation factor pij as extracted from 
the PDF data cr|_ B , er B _ Mg , and from the single atom 

mean-square lattice displacements o~ 2 (i a ) is shown in the 
lower panels of Fig. (open circles) , together with the cor- 
responding correlation factor predicted by the CF model 
(solid lines). In order to extract the experimental value 
of pB-Mg we have taken into account a slight magnesium 
disorder along the c axis, Oq (MgJ = 0.0013 A 2 , m agree- 
ment with the previous analysis of the mean square abso- 
lute displacements a 2 (i a ). We find a positive correlation 
factor for both pb-b ~ 0.1 and ps-Mg ~ 0.1, indicating a 
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slight predominance of the in-phase B-B and B-Mg lattice 
displacements in this experimental probe. Positive values 
of pij are commonly reported in a variety of materials. The 
intuitive explanation is that the in-phase phonon modes 
(acoustic, low optical branch modes) are generally less 
stiff than the opposite-phase optical ones. Our reported 
values of PB-B,PB-Mg ~ 0.1 are however much smaller 
than the correlation factors commonly found in other co- 
valently bonded materials. |35f2D13 1132) This observation 
points out once more that all the motions are decoupled 
and the atoms are behaving largely like independent os- 
cillators. 

Another anomalous feature of MgB2 pointed out by 
this analysis is the lack of a temperature dependence for 
the correlation factors py as compared with the CF model 
and with the standard behavior of other common materials. 
An increase of pij as function of temperature is indeed 
observed in many covalently bonded systems and it is es- 
sentially due to the fact that the thermal population of 
the low frequency in-phase phonon modes is larger than 
the high frequency out-of-phase phonon modes. The lack 
of this temperature dependence in our measurements can 
maybe be attributed to the anharmonic character of the 
high frequency (£?2g) B-B modes and it represents an in- 
teresting anomaly in this material whose physical inter- 
pretation can shed interesting light on the lattice dynam- 
ics in MgB2. Further work on this subject is required. As 
a final point, we can quantify in our model the role of 
the electron-phonon coupling on the correlation factors. 
As shown in Fig. [3] the inclusion of the electron-phonon 
interaction, which leads to a partial softening of the out- 
of-phase Ei g in-plane boron displacements, is reflected in 
a significant reduction Ap^-B ~ —0.03 of the correlation 
factor pb-b while a negligible effect is found on pB-Mg- 

4 Conclusions 

In this work we have investigated the local lattice prop- 
erties of MgB2 paying special attention on the lattice dy- 
namics and the correlations in the B-B and B-Mg first 
neighbor atomic pair motion. We have used the real space 
PDF obtained from high resolution neutron diffraction to 
study the effects of the lattice vibrations on the PDF peak 
widths. The PDF peaks in well ordered crystals such as 
the present case yield important information about the 
underlying atomic potentials through the correlated lo- 
cal lattice dynamics. The data have been modeled using 
both a multi-parameter constant force model and a sim- 
ple Einstein model. We have found that the constant force 
model as well as the Einstein one reproduce the average 
features of the lattice vibrations. This agreement suggests 
that boron and magnesium displacements, both in-plane 
and out-of-plane, are mostly independent of each other. 
The analysis of the PDF peak linewidths permits to evalu- 
ate the correlation for both the nearest neighbor B-B and 
B-Mg atomic pairs. We find a small positive correlation 
factor pb-b ~ 0.1 and ps-Mg ~ 0.1, nearly temperature 
independent, indicating a weakly prevalent in-phase rela- 
tive atomic motion. These results are in contrast with CF 



model which predicts correlation factors increasing with 
the temperature. This discrepancy supports the idea that 
anharmonic effects and strong decay processes for the Ei g 
B bond stretching modes are present, presumably due to 
the strong electron-phonon coupling. 
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